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CD , ABSTRACT 

I We investigate the physical mechanisms that shape the luminosity function of galaxies in hierarchical 

, clustering models. Beginning with the mass function of dark matter halos in the ACDM cosmology, we 

■ show, in incremental steps, how gas cooling, photoionization at high redshift, feedback processes, galaxy 

merging and thermal conduction affect the shape of the luminosity function. We consider three processes 
(f) . whereby supernovae and stellar wind energy can affect the forming galaxy: (1) the reheating of cold disk 

J> ' gas to the halo temperature; (2) expansion of the hot, diffuse halo gas; (3) complete expulsion of cold 

disk gas from the halo. We demonstrate that while feedback of form (1) is able to flatten the faint end of 
the galaxy luminosity function, this process alone does not produce the sharp cut-off observed at large 
luminosities. Feedback of form (2) is also unable to solve the problem at the bright end of the luminosity 
function. The relative paucity of very bright galaxies can only be explained if cooling in massive halos 
is strongly suppressed. This might happen if thermal conduction near the centres of halos was very 
efficient, or if a substantial amount of gas was expelled from halos by process (3) above. Conduction is a 
"■-»« . promising mechanism, but an uncomfortably high efficiency is required to suppress cooling to the desired 

' level. If, instead, superwinds are responsible for the lack of bright galaxies, then the total energy budget 

required to obtain a good match to the galaxy luminosity function greatly exceeds the energy available 
from supernova explosions. The mechanism is only viable if the formation of central supermassive black 
S-^ | holes and associated energy generation play a crucial role in limiting the amount of stars that form in 

the host galaxy. The models that best reproduce the galaxy luminosity function, also give reasonable 
approximations to the Tully-Fisher relation and the galaxy autocorrelation function. 

Subject headings: galaxies: luminosity function — galaxies: formation — conduction — cooling flows 
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& Primack (1999). However, although this form of feed- 
back solves the problem with the faint-end slope of the 
luminosity function (as it was originally designed to do), 
it creates a second problem matching the very sharp (ex- 
ponential with luminosity) cut-off seen at the bright end. 
In addition, the effectiveness with which the gas needs to 
be reheated seems excessive compared to observations of 
galactic outflows (Martin 1999) and to what is seen in sim- 
ulations of this process in individual galaxies (Strickland 
& Stevens 2000). These papers suggest that the mass in 
the outflow is comparable to the gas mass that is turned 
into stars. 

All current models of galaxy formation, calculated us- 
ing either gas-dynamical simulations (Pearce et al. 2001; 
Kay et al. 2002; Katz et al. 2002; Weinberg et al. 2003) 
or semi-analytic techniques (e.g. Kauffmann et al. 1999; 
Somerville & Primack 1999; Cole et al. 2000), exhibit 
strong gas cooling in the central regions of groups and 
clusters. This leads to the formation of extremely bright 
galaxies (which are never seen in reality) unless some addi- 
tional suppression of the gas cooling is assumed. The sup- 
pression mechanisms which have been considered in semi- 
analytical models are: simply switching off cooling and/or 
star formation in the most massive halos, e.g. Kauffmann 
et al. (1999); redistributing the gas within the halo so that 
it becomes so rarefied that its cooling time is longer than 
the age of the Universe, e.g. Cole et al. (2000); or suppress- 
ing cooling until the halo is completely virializcd, e.g. van 
Kampen et al. (1999). The model presented in Cole et al. 
(2000) adopted a low value for the cosmic baryon fraction 
(fib = 0.02), and a model for the halo gas distribution 
in which the core radius was a function of the amount 
of gas that had cooled in previous generations of halos. 
In combination, these assumptions were able to produce 
a good match to the galaxy luminosity function; however, 
both are now disfavored by current observational data (e.g. 
O'Meara et al. 2001; Allen et al. 2001; Ettori et al. 2002). 
If the more recent, higher value of the baryon abundance 
(fib ~ 0.04) is adopted, the growing core radius has little 
effect on the cooling rates (since the cooling radius is then 
significantly beyond the typical core radii). 

In this paper, we investigate physical processes that may 
be responsible for shaping the bright end of the luminosity 
function. In addition to the conventional feedback process, 
in which cold gas in the disk is reheated to the tempera- 
ture of the diffuse gas halo, the feedback energy may be 
used to regulate the formation of the galaxy in two further 
ways. While the wind flowing from the gas disk may con- 
tain a relatively small mass, the energy in the disk outflow 
may be transferred to the existing hot gas halo, causing 
it to expand within the halo potential. This makes the 
central gas more diffuse, lengthens its cooling time and 
so reduces the rate at which gas is supplied to the cold 
disk. This type of process is seen in simulations of the 
effect of injection of relativistic radio-emitting plasma in 
the centres of clusters (Quilis et al. 2001; Bruggen et al. 
2002). An alternative is to assume that the gas expelled 
from the disk does not mix with the virialized diffuse gas 
halo, and that the energy per particle is sufficiently high 
that the material escapes the confining potential and is 
never recaptured. 

A further means by which the supply of cold gas to the 



galaxy disk may be reduced is suggested by recent Chan- 
dra and XMM observations of galaxy clusters (Peterson et 
al. 2001; Tamura et al. 2001; Fabian et al. 2001; Johnstone 
et al. 2002; Bohringer et al. 2002; McNamara et al. 2001; 
Nulsen et al. 2002). These have led to a revival of the 
idea that thermal conduction may be an important source 
of heating in the central regions of clusters (Narayan & 
Medvedev 2001; Gruzinov 2002; Voigt et al. 2002). In- 
deed, the conductivity has been inferred to be close to the 
Spitzer rate expected for a fully ionized plasma (Fabian, 
Voigt & Morris 2002). As first suggested by Fabian, Voigt 
& Morris (2002), heating due to conduction could plausi- 
bly counteract the excessive cooling predicted to occur in 
the most massive halos and thereby explain the dearth of 
highly luminous galaxies in the Universe. 

In the remainder of this paper, we will demonstrate 
quantitatively the importance of each physical mechanism 
and its role in setting the form of the luminosity function. 
We will demonstrate that the milder forms of feedback are 
unable to account simultaneously for the sharp break in 
the luminosity function and for the flat faint-end slope. 
The two processes discussed above, thermal conduction 
and superwinds, may provide the answer to this problem. 
Thermal conduction is capable of suppressing the forma- 
tion of the brightest objects if the conductive efficiency is 
high. Alternatively, a model that includes energetic gas 
expulsion is also able to produce reasonable fits to the 
observed luminosity function, but the energy required to 
power this "superwind" is larger than that available from 
supernova explosions and stellar winds. 

2. THE MODEL 

We compute the luminosity function of galaxies in a cold 
dark matter Universe using a development of the semi- 
analytic model, GALFORM, described by Cole et al. (2000). 
Our changes to the model concern the treatment of the gas 
distribution in halos, and the treatment of feedback. 

2.1. Feedback 

In Cole et al. (2000), the star formation rate was af- 
fected by feedback generated when cold gas is reheated 
by energy injected by supernovae explosions and stellar 
winds. This process transfers gas from the star-forming 
disk of the galaxy to a diffuse corona or halo. The rate at 
which cold gas is reheated is assumed to be related to the 
star formation rate, M+, by: 

Mrcheat- (^) Q M„ (1) 

where Vdisk is the circular speed at the half-mass radius 
of the gas disk, and Vhot and a are adjustable parame- 
ters. Here we include this feedback mechanism, but we 
also consider two further schemes which correspond to 
different processes by which the energy from supernova 
explosions and stellar winds can regulate the gas cooling 
and star formation rates. A further fraction of the total 
energy will be radiated away and so is not available for 
feedback. We assume that the way the energy is appor- 
tioned amongst the three processes is fixed by the local 
supernova environment, and does not depend on the mass 
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of the galaxy or on the instantaneous star formation rate. 1 
Thus, feedback in our model is specified by three param- 
eters: e re heat, fhaio, e S w It is useful to note that the total 
energy available from supernovae is of order 10 49 erg for 
each Solar mass of stars that is formed, for a standard 
IMF. 

1. Disk reheating. This is the feedback scheme ex- 
plored by Cole et al. (2000). Gas is removed from 
the disk of the galaxy at a rate proportional to the 
star formation rate. This gas can be thought of as 
cither being ejected at close to, but below the es- 
cape velocity of the halo, or as being heated to the 
halo virial temperature. In the first case, the gas 
will subsequently be heated back to the halo virial 
temperature, through mixing with the halo gas, or 
through further shocks when a new halo forms. In 
cither case, the gas does not leave the halo and so is 
available for cooling in the future. We parameterize 
the energy invested in reheating the cold disk gas 
as e ro hoat 10 49 erg per Solar mass of stars formed. If 
we approximate the gravitational potential wells of 
galaxies as being self-similar, then the energy re- 
quired to eject the gas from the disk should scale 
as V^ 2 isk , and one obtains formula (1) above with 
a = —2. Specifically, if we assume that gas is 
ejected with energy per unit mass equal to V^ 2 isk , 
then the rate at which gas is reheated is given by: 

(Vdi s k/300kms 1 ) 
We therefore have the following relation between 
Crchoat and Vhot- 

2 



£rehoat — 0.18 



300 km s" 



(3) 



If the circular speed of the disk is equal to Vhot , the 
mass of gas reheated is equal to the mass formed 
into stars. Feedback of this form is motivated by 
the study of Efstathiou (2000), although the effi- 
ciency assumed here is somewhat higher. 

2. Energy injection into the gaseous halo. The feed- 
back energy is assumed to redistribute the hot, dif- 
fuse gas so that the central density of the gas is 
reduced, and the central cooling time is increased. 
The energy injected per unit mass of stars formed 
is parameterized as Chaio 

lO^ergMg 1 . The nature 
of this type of feedback process can be understood 
as follows: the total energy injected into the halo 
grows with time until it becomes comparable to the 
gravitational binding energy of the gaseous halo. 
At this point, the cooling rate drops dramatically. 
Since the binding energy of a halo is proportional to 



halo' the fraction of the total gaseous mass which 
must be turned into stars to achieve this balance 
is strongly dependent on halo mass. Larger halos 
are therefore able to cool a larger fraction of their 
baryon reservoir to form galaxies. This mechanism 
is described in more detail in Bower et al. (2001) 
and, as noted above, is motivated by simulations 



of the injection of radio plasma in cluster centres 
(Quilis et al. 2001; Briiggen et al. 2002). 

3. Gas expulsion from the halo. In this scheme, an 
energy of e sw 10 49 ergMg 1 is invested in heating a 
small mass of the cold gas disk to an energy much 
greater than the binding energy of the halo. If 
the energy is contained within this super-heated 
phase (and is not shared with the diffuse gas halo, 
as assumed in mechanism 2 above), the wind may 
be able to escape completely from the halo. This 
mechanism is motivated by the strong, high velocity 
winds inferred to exist around Lyman-break galax- 
ies (Pettini et al. 2002; Adelberger et al. 2003) and 
has recently been studied theoretically by Shu, Mo 
& Mao (2003). We assume that the energy is in- 
jected with a thermal distribution of particle ener- 
gies, i.e. so that the fraction of mass ejected with 
energy in the range E to E + dE is proportional to 
cxp(— E/E av )dE, where E av is the mean ejection 
energy per particle. The fraction of the supcrwind 
gas with sufficient energy to escape the halo is com- 
puted using this distribution. Some of the expelled 
gas will be recaptured when a deeper potential well 
forms. To estimate this, we calculate the fraction 
of the expelled gas with energy less than the depth 
of the new potential well, and allow this fraction 
to be recaptured. This process continues as deeper 
potential wells form until either all of the gas is re- 
captured, or z = is reached. 

We assume that the mass flux of material in the 
wind as it leaves the disk, M superw ind, is propor- 
tional to the star formation rate, with coefficient of 
proportionality /?sw- The energy in the superwind 
as it leaves the disk is given by 

^supcrwind = 10 49 e sw ^* ergs, (4) 

where 10 49 esw ergs is the energy per unit mass of 
star formation. The characteristic specific energy of 
the wind is simply E m = i?supcrwind/M suporwind = 
jVg W , and so the characteristic velocity of the wind 
is 

V sw = 1002 VW/W km/s- (5) 

We assume that material ejected from the disk re- 
quires an energy per unit mass AswVj isk in order 
to escape also from the halo. Therefore, in the 
case in which all of the wind material has the same 
energy, superwinds will be driven from halos with 
Vdisk < 709 y/ esw / /?sw Asw km/s, and not from 
deeper potential wells. For typical halos in our 
model Asw = 2.9 gives a good estimate of the en- 
ergy required to escape the combined gravitational 
pull of the galaxy and its dark matter halo, and 
so we use this value throughout this work. If we 
assume a distribution of energies in the wind, then 
there will be a smooth transition, as Vdisk increases, 
between the regime in which the wind escapes and 
the regime in which it is retained. Thus, the actual 



1 In principle the supernova's local environmcnt- 
but we neglect any such dependence here. 



-the interstellar medium — could depend on global properties of the galaxy, such as its mass, 
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mass flux in the supcrwind escaping from the halo 
is given by 

Msupcrwind = /eject (^W AswVd| sk ) /3swAf*, (6) 

where, for our chosen energy distribution, 
/eject (^) — cxp(— x). The material which does 
escape may be recaptured by larger halos forming 
at later times. Our estimate of the fraction of the 
wind which escapes and that which can be recap- 
tured by larger halos is described in Appendix A. 

We regard this mechanism as being more uncertain 
than mechanisms 1 and 2 because the expelled gas 
must leave the halo without sharing its energy with 
the diffuse hot component. It must therefore punch 
a well collimated hole through the halo, or remain 
contained within a buoyant bubble which is con- 
verted out of the halo (Quilis et al. 2001; Springcl 
& Hernquist 2003; Kay et al. 2002). 

In general, it seems plausible that some fraction of the 
total feedback energy will be processed into each one of 
these forms, and we consider models which include a com- 
bination of these processes on an equal footing with models 
that include only one. 

2.2. Halo Gas Distribution 

The hot gaseous component in dark matter halos is as- 
sumed initially to have a density profile, p g (r), given by 
the /3-model, ie, at radius r in the halo (whose dark mat- 
ter density profile is assumed to have the NFW profile, 
Navarro, Frcnk & White 1996, 1997), 

Ps{r) = [l + (r/ P r c ) 2 p/2' (?) 

where po is the density at the centre of the halo, r c is the 
radius of the "core" and (3 is a parameter which sets the 
slope of the profile on scales larger than r c . 

Departing from the prescription of Cole et al. (2000), 
we adopt a gas density profile in the absence of energy in- 
jection with fixed r c = 0.07r v ; r (and (3 — 2/3) in all halos. 
This provides a reasonable match to gas-dynamic simula- 
tions of non-radiative gas in clusters (e.g. Eke, Navarro 
& Frenk 1998) and to the observed X-ray profiles of re- 
laxed clusters (e.g. Allen et al. 2001). The normalization 
of the profile is determined by the total diffuse gas mass 
remaining in the halo, and the temperature of the gas is 
set assuming hydrostatic equilibrium. As a boundary con- 
dition, we set the temperature at the virial radius equal 
to the virial temperature (eqn. 4.1 in Cole et al. 2000). 
This "default" profile is modified if the diffuse gas gains 
further energy ( "excess energy" ) as a result of energy in- 
jection (process 2 in Section 2.1). We use the prescription 
described in Bower et al. (2001), in which the excess en- 
ergy first causes the slope of the gas profile to decrease 
down to a minimum value of /3 m j n = 0.2, after which it 
increases the boundary temperature of the gas halo. The 
effect of the excess energy is to decrease the central density 
of the gas, lengthening its cooling time. Mass is conserved 
by pushing some of the diffuse gas outside the halo; how- 
ever, this gas can be recaptured as the total halo mass 
(and thus the gravitational binding energy) increases. We 
assume that the excess energy is conserved during mergers 
between halos, although the results are not qualitatively 



affected by a small dilution or amplification of energy dur- 
ing halo mergers. 

We define the effective cooling time at radius r, t' cool (r), 
as the maximum true cooling time (i.e. that defined by 
Cole et al. 2000) occurring at smaller radii, plus the free- 
fall time from r to the halo centre. This ensures that 
dt' cool /dr > 0, so the cooling radius is always a smoothly 
increasing function of radius. Experiments with different 
approaches show that the results we present here are not 
sensitive to the details of this prescription. 

2.3. Conduction 

Conduction in the ionized gas can transport energy into 
the inner regions of the halo, effectively increasing the cool- 
ing time of the gas there. The rate at which energy is de- 
posited into the shell between radii r and r + dr is given 
by: 

Sdr = 47r-i ( nr 2 ^-) dr, (8) 
dr \ dr J 

where the conductivity, k, may depend on radius through 
its temperature dependence. We approximate £ as 

£ = a C ond47TK s T, (9) 

where k s is the Spitzer conductivity (Spitzer 1962), and 
a con d is a parameter that absorbs the dependence on the 
shape of the temperature profile as well as any difference 
between the actual conductivity and the Spitzer rate. For 
a power-law temperature profile, T a r a , and conductiv- 
ity of the Spitzer form, n s cx T 5 / 2 , a con d = fs P o(l + 7a/2) 
where / sp is the ratio of the conductivity to the Spitzer 
value. Adopting a temperature profile with a = 0.4 as 
suggested by recent X-ray observations (Voigt et al. 2002) 
gives a C ond = 0.96/ S p. Adopting a linear temperature gra- 
dient gives a co „d = 4.5/ sp . 

The heating rate due to conduction is subtracted from 
the radiative cooling rate to give a net cooling rate for the 
gas. This net cooling rate is used to compute the cool- 
ing time. Conduction causes the cooling radius to become 
smaller than in the standard model. The result is a sup- 
pression of cooling in hot halos, and the mass at which this 
effect becomes important is determined by the parameter 

^cond- 

2.4. Photoionization and Merging 

While we employ the detailed model of galaxy merg- 
ing developed by Benson et al. (2002), we choose to use a 
much simpler model of the effects of reionization than used 
in that paper, in order to incorporate them easily into our 
calculation of the galaxy luminosity function. We simply 
assume that galaxy formation is completely suppressed by 
reionization in dark matter halos with circular velocities 
below Vroion after rcdshift z ro j on . Unless otherwise stated, 
we adopt V Ic i on = 50kms _1 and 2; ro i on = 6. With this 
choice, this simple model matches the results of the full 
calculation of Benson et al. (2002) quite well. 

3. RESULTS 

Throughout this paper, we use the same parameter val- 
ues adopted by Benson et al. (2002), with the exception 
of a larger baryon fraction corresponding to = 0.045, 
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consistent with constraints from Big Bang nucleosynthe- 
sis (O'Meara et al. 2001). We assume a ACDM universe 
with mean matter density f^o = 0.3, cosmological constant 
term fi A = 0.7, Hubble constant 2 H = 70km s _1 Mpc -1 
and linear fluctuation amplitude on spheres of radius 8h 
Mpc, cr 8 = 0.9. Of these parameters, the uncertainty in the 
value of er 8 has most impact on the model results. While 
several studies support a value of as ~ 0.9 (e.g., Bacon et 
al. (2002), Hoekstra et al. (2002) [gravitational lensing]; 
Eke, Cole & Frenk (1996), Vianna et al. (2002) [cluster 
abundance]; Spergel et al. (2003) [cosmic microwave back- 
ground (CMB)]; Sievers et al. (2003) [Sunyaev-Zeldovich 
effect]), other recent analyses have suggested lower val- 
ues, as ~ 0.7 (e.g., Peacock (2003) [large scale structure]; 
Melchiori et al. (2002) [CMB]; Jarvis et al. (2003) [grav- 
itational lensing]; Allen et al. (2003), Smith et al. (2003) 
[cluster abundance]). A discussion of recent results may 
be found in Wang et al. (2003). Except where specified, we 
show models based on er 8 = 0.9; however, (as we shall see) 
taking the lower value, <r 8 = 0.7, considerably eases the 
energy budget and/or reduces the conduction efficiency re- 
quired to match the galaxy luminosity function. Through- 
out, we use the halo mass function derived from N-body 
simulations by Jenkins et al. (2001) 3 , instead of the Press- 
Schcchtcr mass function used by Cole et al. (2000). 



odel 1 = Halo mass function 
Model 2 = Cooling 
Model 3 - Photoionization 
Model 4 = Merging 




Fig. 1. — The K-band luminosity function of galaxies. The points 
show the observational determinations of Cole et al. (2001, circles), 
Kochanek et al. (2001, squares) and Huang et al. (2002, z < 0.1, 
stars). Lines show model results. Model 1 (dashed line) shows the 
result of converting the dark matter halo mass function into a galaxy 
luminosity function by assuming a fixed mass-to-light ratio chosen 
to match the knee of the luminosity function. Model 2 (dotted line) 
shows the result from GALFORM when no feedback, photoionization 
suppression, galaxy merging or conduction are included. Models 3 
and 4 (long dashed and solid lines respectively) show the effects of 
adding photoionization and then galaxy merging. 

We compare our model with recent determinations of 

2 Here and below h denotes the Hubble constant in units of 100km s _1 Mpc -1 

3 Note that we use the Jenkins et al. (2001) mass function to compute the abundances of halos at z = 0, but retain the Press-Schechter 
approximation when computing halo merger trees. An improved calculation would use a self-consistent computation of halo merging histories. 

4 We set X = 1 in the notation of Cole et al. (2000) 



the K-band luminosity function assuming a Kennicutt stel- 
lar initial mass function (Kennicutt 1983). In order to 
facilitate comparison between models, we have kept the 
IMF fixed, and assumed a negligible fraction of brown 
dwarf stars 4 . We choose the K-band in order to mini- 
mize the sensitivity of our results to recent star formation 
and to dust obscuration. The model of Cole et al. (2000) 
includes a detailed and fully self-consistent calculation of 
dust extinction which is used in this work. We find, how- 
ever, that dust-obscuration has a negligible effect on our 
results for the K-band luminosity function (typically shift- 
ing the bright end faintwards by less than one tenth of a 
magnitude.) For the observational comparison, we use the 
local K-band luminosity functions measured by Cole et 
al. (2001) and Kochanek et al. (2001) (both based on the 
2MASS survey) and the local (z < 0.1) luminosity function 
derived from the much deeper K-band survey of Huang et 
al. (2002). The analysis by Huang et al. suggests a faint- 
end slope («k = —1-37 ± 0.10), steeper than the values 
found by Cole et al (ax = —0.93) and by Kochaneck et al. 
(ax = —1-09). These latter two are also in good agreement 
with the faint end slope of the z-band luminosity function 
measured by Blanton et al. (a z = —1.08) from the SDSS 
survey. The z-band data should also be little affected by 
residual star formation and dust extinction, but have a 
deeper surface brightness limit. These discrepancies indi- 
cate that there remain significant systematic uncertainties 
in current measurements of the faint end of the K-band 
luminosity function, perhaps due to luminosity from low- 
surface brightness regions of galaxies being missed as re- 
cently suggested by Andreon (2002). 

3.1. A Model with constant mass-to-light ratio 

In Fig. 1 we show the simplest possible model of the 
luminosity function which we call Model 1 (shown as the 
dashed line). In this model, the mass function of dark 
matter halos (Jenkins et al. 2001) has been converted into 
a luminosity function simply by assuming a fixed mass-to- 
light ratio (M/Lk — 11M q /Lk,©), chosen so as to match 
the knee of the observed luminosity function. As is we'll 
known, this produces a luminosity function which is much 
steeper at the faint end than is observed, and also fails to 
cut off at bright magnitudes (the halo mass function does 
possess a cut-off, but it occurs at much higher mass and 
lower abundance than shown in the plot). 

3.2. Model including Cooling only 

White & Rees (1978) argued that the difference between 
the halo mass and the galaxy luminosity functions is due 
to the dependence of the gas cooling time on halo mass and 
to feedback processes. We use the galform semi-analytic 
model to follow gas cooling and star formation in a merging 
hierarchy of dark matter halos in the ACDM cosmology. In 
order to illustrate the simplest possible model first, we do 
not include photoionization suppression, feedback, galaxy 
merging or conduction. The result is Model 2 (shown as 
a dotted line in Fig. 1). It clearly displays the "overcool- 
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ing problem": gas has cooled into the smallest halos re- 
solved in the calculation, producing an overabundance of 
faint galaxies. As a result, the faint-end slope of the lu- 
minosity function is much too steep. It should be noted 
that the results for this cooling-only model are sensitive 
to the mass resolution of the GALFORM merger trees. The 
cooling time for halo gas decreases with decreasing halo 
mass down to halo virial temperatures around 10 4 K. Be- 
low this scale, cooling becomes ineffective unless molecular 
hydrogen is abundant. This means that for the luminosity 
function to be fully converged at all galaxy luminosities, 
we would need to resolve halos with virial temperatures as 
low as 10 4 K at all redshifts at which there is significant 
cooling. To obtain the results shown for Model 2, we ran 
GALFORM at the highest mass resolution that was com- 
putationally feasible, which comes close to resolving 10 4 K 
halos at z = 0. We believe that our results are substan- 
tially converged, but we cannot be certain that they would 
not change if the mass resolution were increased further. 
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Fig. 2. — Model 5. Starting from Model 4, disk-reheating is 
added in order to suppress the formation of small galaxies. Results 
are shown for three levels of energy input (^reheat = 0.03, 0.13 and 
0.41). The data points are the same as Fig. 1. 



3.3. Model with photoionization and merging 

The problem of mass resolution encountered in Model 2 
is effectively removed once we include feedback processes 
which are effective in low-mass halos. Model 3 (long- 
dashed line in Fig. 1) shows the effects of including pho- 
toionization suppression. As described by Benson et al. 
(2002), the formation of low-mass galaxies is suppressed 
after the Universe is reionized and, as we noted in §2, we 
adopt a simplified model of this process which nevertheless 
is a good approximation to a full calculation. The resolu- 
tion of the GALFORM merger trees is then sufficient to fol- 
low the lowest mass halos that are able to form galaxies, re- 
sulting in a converged solution. However, photoionization 
suppression alone is unable to produce a sufficiently flat 
faint end in the luminosity function. In addition, Model 3 
also produces more bright galaxies, since less gas has been 
locked into stars in the smallest halos. When both galaxy 
merging (which was artificially switched off in Models 1 
and 2) and photoionization are included, as in Model 4 
(solid line in Fig. 1), the faint end remains too steep and 
even more bright galaxies are produced. In the models 
that follow, we investigate the effects of including feedback 
processes in addition to photoionization and merging. 



3.4. Model with feedback — disk reheating 

A solution to the faint-end problem is illustrated by 
Model 5.3 in Fig. 2. Here, feedback is included through 
the reheating of disk gas in star-forming galaxies. We 
use the standard prescription of Cole et al. (2000), but 
with a larger value of e ro hoat = 0.41 (equivalent to Vhot = 
450 km s -1 ) which is required in order to obtain a similar 
faint end slope to that in Cole et al. for the larger value of 
fib assumed in this work. This form of feedback flattens 
the luminosity function considerably, resulting in reason- 
ably good agreement with the observed faint end. While 
the slope is not as flat as that measured by Cole et al. 
(2001) or Blanton et al. (2003), it is in good agreement 
with the steeper slope reported by Huang et al. (2002). 
This achievement carries a price, however — the overabun- 
dance of bright galaxies (formed through excessive cooling 
in massive halos) is exacerbated, as there is now a much 
greater mass of diffuse hot gas remaining in the larger 
halos. This gas is sufficiently dense that the central re- 
gions are able to cool; consequently, Model 5.3 produces 
far too many bright galaxies. This result depends little on 
the choice of er 8 . Adopting as = 0.7 makes the brightest 
galaxies only 0.5 mag fainter. This clearly demonstrates 
a longstanding problem in semi-analytic models: previous 
calculations have either assumed low values of f2b, or have 
invoked rather artificial ways to prevent the cooling which 
forms these over luminous objects. 
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gas from dark matter halos at such high energies that it is 
subsequently unable to cool. 
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Fig. 3. — Model 6. These models illustrate the effect of energy 
injection. Starting from Model 5.3, the effect of heating the diffuse 
gas halo is included. Models 6.1-6.3 have ehalo = 0.1, 0.2 and 1.0 
respectively, in addition to e r chcat = 0.41. Increasing the energy 
available for this form of feedback suppresses the formation of both 
bright and faint galaxies, but does not produce a sharp break in the 
luminosity function. 

The energy requirements of Model 5.3 are substan- 
tial, but not excessive. The reheating energy of 0.41 x 
10 49 erg Mq 1 should be compared to the total energy avail- 
able from supernovae explosions which is approximately 
0.7x 10 49 erg Mq 1 for a Salpeter IMF, or 0.9 x 10 49 erg Mq 1 
for the Kennicutt IMF adopted in GALFORM. In a halo of 
circular velocity 250 km/s, the mass of gas reheated is 
more than 3 times the mass of gas formed into stars. If 
the level of reheating is reduced, as illustrated by mod- 
els 5.1 (e re heat = 0.03) and 5.2 (e rohcat = 0.13), then the 
formation of small galaxies is not suppressed sufficiently 
to match the observed luminosity function. 

3.5. Models with energy injection 

In Models 6.1-6.3 (Fig 3), we investigate the effect of 
assuming that a fraction of the supernova and stellar wind 
energy heats the diffuse gas halo, causing it to expand. 
This form of feedback suppresses the formation of both 
bright and faint galaxies, but it does not produce a sharp 
break in the luminosity function. Models 6.1 and 6.2 illus- 
trate how as the energy spent in heating the diffuse halo 
increases, the break in the luminosity function becomes 
less pronounced. If the heating is made even stronger, 
the resulting luminosity function approaches a power-law 
(Model 6.3). This form of feedback clearly cannot solve 
the problem of overproduction of bright galaxies. 

3.6. Model with conduction 

We now consider two possible schemes that are capa- 
ble of producing a good match to the observed luminosity 
function. The first involves balancing radiative cooling 
with thermal conduction. The second involves expelling 
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Fig. 4. — Model 7. These models illustrate the effect of thermal 
conduction. In Model 7.3 (a con d = 25), conduction is assumed to 
be highly efficient (it is unlikely that such a high efficiency is physi- 
cally plausible). More realistic conduction efficiencies arc illustrated 
in Models 7.2 (o cond = 1) and 7.1 (a con d = 0.1). For Model 7.4, 
we adopt a lower value for as; a conduction efficiency of Cf cont j = 7 
then gives a reasonable match to the observed luminosity function. 
In all cases, the energy feedback parameters parameters are set to 
^halo = 0.0 and e rchoat = 0.41. 



Thermal conduction is expected to imprint a special 
scale on the galaxy population because of the strong 
temperature dependence of the Spitzer conductivity rate. 
Model 7.3 (solid line in Fig. 4) shows the result of including 
conduction with a con d = 25 in Model 5.3 (e ro h at = 0.41, 
£haio = 0) . Such a high value of the conductivity is indeed 
effective in suppressing the formation of the most massive 
galaxies since it prevents efficient gas cooling in group and 
cluster-sized halos. The brightest galaxies in the luminos- 
ity function are instead built through mergers. The result 
is a rather good match to the observed galaxy luminosity 
function. However, the conduction efficiency that we have 
assumed is extremely high. To operate at the required 
level, we must assume both that the conductivity is not 
suppressed below the Spitzer value, and that the effective 
temperature gradient is steeper than T oc r 2 5 in the re- 
gion of the cooling radius. Note that the Spitzer formula 
for thermal conductivity in an ionized plasma breaks down 
if the conductivity becomes too high, i.e. if the conduc- 
tion "saturates" as described by Cowie & McKee (1977). 
The models shown in this paper do not take account of 
this saturation limit. However, using the estimate of the 
saturated heat flux from Cowie & McKee (1977), we have 
checked that our results are not significantly affected by 
saturation (for model 7.3, which has the most extreme con- 
duction of all our models, there is only a small increase in 
the number of the very brightest galaxies). 

In Models 7.1 and 7.2, we show the effect of assuming 
a more modest conduction efficiency (a CO nd = 1.0 and 0.1 
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respectively). In these models, conduction is not sufficient 
to suppress cooling in the larger halos adequately. 

If we adopt a lower value for as, however, a lower con- 
duction efficiency gives a reasonable match to the observed 
luminosity function. Model 7.4 shows the luminosity func- 
tion for the case as = 0.7 and a con d = 7. This conduction 
efficiency could be achieved if the temperature gradient 
was T oc r 1 - 3 and the conduction was only slightly sup- 
pressed below the Spitzer value. Although this is still a 
high rate of conduction, it offers a promising route for ex- 
plaining the bright end of the galaxy luminosity function. 
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Fig. 5. — Model 8. These models illustrate the effect of super- 
winds. In Model 8.1, an energy of e sw = 0.27 drives a weak super- 
wind (with /3gw — 3); disk reheating has efficiency e re h eat = 0.13 
and there is no heating of the diffuse halo (ehalo = 0)- A much 
more powerful wind is needed to create a break in the luminosity 
function. Model 8.2 (e sw = 5.0; /3sw = 1) illustrates the effect of 
increasing the superwind power. An improved match to the lumi- 
nosity function can be achieved with the same power if ag is lower. 
Model 8.3 shows a model with e sw = 5.0, /3sw = 1 an d °"8 = 0.7. 



3.7. Model with superwinds 

The expulsion of gas from halos at high energy can, in 
principle, strongly suppress the formation of later gener- 
ations of galaxies, hence affecting the shape of the lumi- 
nosity function. Starting from Model 5.2, we add further 
feedback energy that expels cold gas completely, not only 
from the disk but also from the halo. The superwind must 
have high energy in order that the expelled material not 
be recaptured by more massive halos. The effect of a low 
power superwind is illustrated by Model 8.1 (dashed line 
in Fig. 5). This model, with e sw = 0.27 and /?sw = 3, has 
a relatively weak superwind. This gas expulsion is in ad- 
dition to the reheating of cold disk gas (e r ohcat = 0.13); we 
have assumed that there is no heating of the diffuse halo 
(ehaio = 0.0). Although the winds eject a large amount 
of gas, most of the material is recaptured as larger ha- 
los collapse and the luminosity function differs little from 
Model 5.2. 

In Model 8.2 (dotted line), we have set e sw = 5.0 and 



Psw = 1, corresponding to a mean energy per superwind 
particle of E av — 15keV. Such an energetic wind is re- 
quired to ensure that very little material is recaptured by 
group halos. In this model, the superwind dominates the 
feedback energy budget; indeed, the total energy required 
(5.13 x 10 49 ergMQ 1 ) significantly exceeds that available 
from supernovae alone. The model comes much closer to 
matching the luminosity function, but still overproduces 
bright galaxies. We can increase the superwind energy 
further, but a factor of 2 increase only results in a small 
improvement in the luminosity function. If we include 
conduction as well as superwinds, it is, of course, possi- 
ble to improve the match but a high conduction efficiency 
(cKcond ^ 1) is still needed. Increasing the mass loading of 
the wind substantially (by increasing /?sw) results in too 
few galaxies around the knee of the luminosity function. 

As we found in the case of conduction, the luminosity 
function can be matched more easily if we adopt a lower 
value for as- The case e sw = 5.0, /?sw — 1, ^8 = 0.7 is 
illustrated by Model 8.3 (solid line). Given the uncertain- 
ties of our recapture prescription, this model gives a rea- 
sonable match to the luminosity function; it has a strong 
break at the correct luminosity, and overproduces bright 
objects only marginally. There are a variety of ways to 
further improve the match to observations: we could in- 
crease the energy injection (but e sw > 10.0 is required) or 
increase the mass loading of the wind so that the curve is 
shifted faintwards. An alternative strategy is to allow for a 
low level of conduction: a con d ~ 1 is sufficient to produce 
a significant improvement in the match to the luminosity 
function when as = 0.7 and superwinds are present. 

The parameters a con d and egw are highly degenerate 
in their effects on the luminosity function (i.e. increasing 
either suppresses the bright end). While current computa- 
tional limitations make it prohibitively expensive to per- 
form an accurate % 2 fit of the model parameters to the 
data, a crude estimate of the \ 2 surface for these two pa- 
rameters shows that the data prefer models with strong 
superwinds (esw ~ 6) and high conductivity (a C ond ~ 30) 
for a model with as — 0.93. Lowering a s to 0.7 reduces the 
requirements to esw ~ 3 and a con d w 25. Further inves- 
tigation of the x 2 surface would require a Bayesian prior 
to specify formally a physically allowed range for this pa- 
rameter. 

3.8. Other Considerations 

While our primary goal in this paper is to examine 
the luminosity function of galaxies, it is prudent to check 
whether our models are in reasonable agreement with 
other basic properties of the galaxy population such as 
the Tully-Fisher relation and the galaxy autocorrelation 
function. We compare the models that best fit the galaxy 
luminosity function to these observables. We retain the 
model parameters found earlier and we do not attempt to 
adjust these or any other parameters in this comparison. 
We defer a more detailed comparison of our models with a 
wider range of observational constraints to a future paper. 

3.8.1. Tully-Fisher Relation 

Simultaneously matching the galaxy luminosity func- 
tion and the Tully-Fisher relation has been a long-standing 
problem for CDM-based semi-analytic models (see, for ex- 
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ample, White & Frcnk 1991, Kauffmann, White & Guider- 
doni 1993 and Cole et al. 2000). In particular, Cole et al. 
(2000) found that while their best model reproduced the 
observed slope and scatter in the Tully-Fisher relation, the 
predicted circular velocities for galaxy disks were about 
30% larger than the values measured for the data. 

Figure 6 shows the Tully-Fisher relation predicted by 
four of our models (including those which most successfully 
match the observed luminosity function). It is readily ap- 
parent that our models all perform reasonably well in that 
they reproduce the slope and scatter of the observed rela- 
tion, although they do not match the zero-point exactly. 
Model 8.3 does best, being offset by only approximately 
10% in circular velocity over most of the range shown. 
This model, in fact, performs rather better than that of 
Cole et al. (2000), although an even better match would 
be desirable. This may involve revising the assumption of 
adiabatic invariance made in calculating the effects of the 
disk on the inner structure of the dark matter halo (see 
Cole et al. 2000) or considering the effects of the angular 
momentum of the superwind material or the details of our 
star formation prescription. 
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Fig. 6. — The I-band Tully-Fisher relation of galaxies. Circles 
show the observational sample of Mathcwson, Ford & Buchhorn 
(1992). Lines show the median relation for model galaxies. We 
follow Cole et al. (2000) and select only galaxies with a bulgc-to- 
total ratio in dust-extinguished I-band light between 0.02 and 0.24 
and with a disk gas fraction of at least 10%. Circular velocities for 
model galaxies are computed at the disk half-mass radius. Errorbars 
enclose the inner 80% of the distribution of model galaxies. 

3.8.2. Galaxy Correlation Function 

One of the most remarkable successes of semi-analytic 
models of galaxy formation in the ACDM cosmology is 
their ability to match the observed galaxy autocorrela- 
tion function (Kauffmann et al. 1999; Benson et al. 2000). 
As pointed out by Benson et al. (2000), the clustering of 
galaxies may be understood physically by reference to the 
halo occupation distribution (see, for example, Peacock 
& Smith 2000, Cooray & Sheth 2002), which specifies 
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the number of galaxies that occupy dark matter halos of 
a given mass. We show, in Fig. 7, the mean number of 
galaxies brighter than Mb — 5\ogh = —19.5 per halo as a 
function of halo mass for four of our models. For reference 
we also plot the fits obtained by Bcrlind et al. (2003) to the 
semi-anaytic model of Cole et al. (2000) and to a smooth- 
particle hydrodynamics calculation of galaxy formation as 
solid and dotted lines respectively. The two models con- 
sidered by Berlind et al. (2003) both gave reasonable fits 
to the galaxy correlation function. 
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Fig. 7. — Symbols show the mean number of galaxies brighter 
than Mb — 51og/i = —19.5 as a function of halo mass for some of 
our models. For reference, lines show the fitting functions of Berlind 
et al. (2003). The solid line shows their fit to the results of the Cole 
et al. (2000) semi-analytic model while the dotted line shows their 
fit to a smoothcd-particlc hydrodynamics calculation. 



Our new models give rise to a mean halo occupancy 
which is qualitatively similar to those found by Berlind 
et al. (2003) - a power-law at high masses, which flattens 
before cutting off sharply at some critical mass. Quantita- 
tively, the models with conduction (7.3 and 7.4) produce 
significantly fewer galaxies in low and high mass halos, 
while agreeing with the Berlind et al. (2003) fits at a halo 
mass of a few times lO 12 ft. -1 M0. The models with super- 
winds (8.1 and 8.3) agree more closely with the Berlind 
et al. (2003) fits, although model 8.3 produces somewhat 
fewer galaxies in clusters. Since none of our models dif- 
fer greatly from those of Bcrlind et al. (2003), we may 
expect them also to produce correlation functions in rea- 
sonable agreement with the observational data, at least 
on intermediate and large scales. On smaller scales, the 
mean number of galaxies per halo is not sufficient to spec- 
ify the correlation function; this requires knowledge of the 
distribution of galaxy pairs per halo. 
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Fig. 8. — The galaxy real-space two-point correlation function. 
Circles show the observational determination of Baugh (1996). The 
heavy solid line shows the correlation function of dark matter for 
our adopted ACDM cosmology (with erg = 0.9). The remaining lines 
show results for model galaxies brighter than M-g — 5 log h = —19.5, 
as indicated in the legend. 

Figure 8 shows correlation functions (computed using 
the techniques of Benson et al. 2000) for the same four 
models, compared to both the predicted dark matter cor- 
relation function and the observed galaxy correlation func- 
tion. As expected, all models do well at matching the 
data on scales larger than about lft, _1 Mpc. On smaller 
scales, all models are 'antibiassed' with respect to the 
dark matter, but they produce somewhat too much clus- 
tering. Benson et al. (2000) found similar results for their 
semi-analytic model. Model 8.3 gives, in fact, very simi- 
lar clustering results to theirs (for the same ACDM cos- 
mology). The other models produce even stronger correla- 
tions on small scales, providing a worse match to the data. 
Model 8.3 does better than Model 8.2 because it places 
relatively fewer galaxies in clusters, thereby having many 
fewer pairs on small scales. We speculate that a model 
which performs even better at matching the observed lumi- 
nosity function (e.g. producing a sharper cut-off at bright 
magnitudes) may also do even better at matching the ob- 
served correlation function. 

4. DISCUSSION 

We have examined the key physics thought to be neces- 
sary to explain the shape of the galaxy luminosity function 
in a cold dark matter universe. While the cooling and con- 
densation of gas in a merging hierarchy of dark matter ha- 
los remains the fundamental process through which galax- 
ies form, we have demonstrated that at least two other 
processes must act to shape the luminosity function. A 
minimal model requires the inclusion of feedback mecha- 
nism^) (beyond the heating resulting from the photoion- 
ization of the pregalactic gas) to flatten the faint end of the 
luminosity function and to suppress cooling at the centres 
of the massive halos of groups and clusters. If a fraction of 



the energy liberated by supernovae and winds goes into re- 
heating disk gas and/or heating the diffuse gas halo, then 
the formation of faint galaxies is suppressed, resulting in 
a flattened faint-end slope that matches the available ob- 
servational data adequately. These mechanisms on their 
own, however, are unable to produce a sharp cut-off at 
the bright end of the luminosity function. We have shown 
that there are two possible (but quite extreme) processes 
that can achieve this: (1) thermal conduction at about 
or above the Spitzer rate; (2) expulsion of gas from halos 
in superwinds at temperatures high enough to prevent its 
subsequent recapture. 

The high value of a con d required to suppress the bright 
end of the luminosity function is discouraging since it 
implies both that the conductivity must be close to the 
Spitzer value (despite the presence of \iG magnetic fields 
in clusters; Narayan & Medvedev 2001; Taylor et al. 2002), 
and that the effective temperature gradient must be some- 
what steeper than the gradients observed in galaxy clusters 
(Allen et al. 2001; Fabian, Voigt & Morris 2002)'. 

However, our method for calculating the effect of con- 
duction is highly simplified. In particular, there are two 
issues that are not well addressed. Firstly, our calcula- 
tion is based on the heat flux through a shell. We have 
not considered the total extent of the region in which 
conduction must be effective. We can define a radius 
'"out, such that the initial thermal energy in the region 
'"cool < r < r out is equal to the total energy radiated from 
the region r < r coo i, where r coo i is the cooling radius in 
the presence of conduction. Thermal energy needs to be 
conducted from a radius at least as large as r out in or- 
der to balance the radiative losses from smaller radii. For 
an isothermal halo with an r~ 2 density profile, we find 
^cooi^out = r'^ ool , where r' cool is the cooling radius calcu- 
lated in the absence of conduction. The radius r out can 
be large, but the temperature gradient that is required to 
transport heat effectively across this region is a strongly 
declining function of radius. Thus, so long as r out lies 
within the virialized region of the cluster, this is unlikely 
to modify our solution significantly. Secondly, our method 
assumes that the logarithmic temperature gradient within 
'"cool is constant. If, instead, the gradient is peaked in 
this region, our method would underestimate the heat de- 
posited in the shell at r coo i. These two simplifications 
(which affect the required conduction coefficient in oppo- 
site senses) are difficult to model accurately and without 
ad-hoc assumptions. Clearly, this is a problem that needs 
to be addressed by numerical simulations of cooling in con- 
ductive plasmas. These would allow us to calibrate our 
simple model of conductive heating and to infer what con- 
ductivity is required, by removing the degeneracy between 
the shape of the halo temperature profile and the suppres- 
sion factor f sp for the conductivity. 

The alternative model invokes highly energetic "super- 
wind" outflows. The expulsion of gas from the halo reduces 
the reservoir of baryons available for cooling at later times, 
and reduces the number of stars formed by each parcel of 
cold gas. However, the fraction of baryons required to be 
in this "super-heated" phase at z = is quite high, and 
the energy budget is exceptional. We found the energy 
budget of a promising superwind model to be in region of 
5 x 10 49 ergs per Solar mass of stars formed. It is quite 
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possible our model may not treat the recapture process 
adequately and that we may therefore be overestimating 
the total energy requirement. It is difficult to model the 
process better, however, because the expelled gas may cool 
adiabatically and then be accreted by a completely differ- 
ent halo; alternatively, some material may escape into void 
regions and evade recapture altogether. A further prob- 
lem is the need to assume a particular distribution for the 
wind particle energies. Our assumption that the distri- 
bution is thermal may lead us to overestimate the energy 
requirement. 

On average, the mass of baryons ejected from halos 
amounts to 7% of mean cosmic value. This corresponds to 
about 1.3 times the fraction of baryons turned into stars. 
A compilation of observations by Martin (1999) suggests 
that the total mass loss rate from disk galaxies could cer- 
tainly be as high as twice the star formation rate. If all of 
this mass loss were energetic enough to produce a super- 
wind, these observations would provide some justification 
for the high value of e sw required in our model. The total 
energy budget, however, exceeds that available from su- 
pernovae (~ 0.9 x 10 49 ergM~ 1 for the Kennicutt IMF) 
even though we have not allowed for radiative losses. Al- 
though the exact energy output from supernovae is uncer- 
tain, and the conditions we have assumed for superwind 
escape are also uncertain, this suggests that the supernova 
energy must need to be supplemented by energy released 
during black hole formation (Enfilin et al. 1998; Bower et 
al. 2001; Cavaliere et al. 2002). 

Finally, we have briefly examined the ability of some of 
our most successful models to match the observed Tully- 
Fishcr relation and galaxy correlation function. With- 
out adjusting any parameters, we find, as in many previ- 
ous semi-analytic calculations, that our models adequately 
match the slope and scatter of the Tully-Fisher relation, 
but tend to overpredict the circular velocity at a given lu- 
minosity (or, equivalently, underpredict the luminosity at 
given circular velocity.) In most cases, the agreement is 
at least as good as in previous studies, but for one of our 
superwind models the discrepancy in the circular velocity 
zero-point is reduced to only 10%. At least one of our mod- 



els (Model 8.3) matches the galaxy autocorrelation func- 
tion reasonably well (as well as the semi-analytic model 
of Benson et al. (2000)). Importantly, this success reflects 
the fact that clusters in this model contain relatively fewer 
galaxies, and thus many fewer pairs on small scales, than 
clusters in other models. This demonstrates how mea- 
surements of galaxy clustering place important constraints 
on halo mass-dependent feedback processes such as those 
considered in this study. Our results encourage a more de- 
tailed investigation of how our models fare when compared 
with a wider range of observables, a programme that we 
intend to pursue in future work. 

5. CONCLUSIONS 

One of the simplest properties of the galaxy population, 
the luminosity function, is now reasonably well determined 
observationally (although, as we have noted, some sig- 
nificant disagreements about the faint-end slope remain). 
We have shown that our theoretical understanding of this 
fundamental property has progressed enormously in re- 
cent years to the extent that it is now possible to explain 
the form of the luminosity function in quantitative de- 
tail. Nevertheless, two crucial ingredients — feedback and 
conduction — are understood only poorly at best. Further 
study of these physical processes is required before we can 
truly claim a satisfactory understanding of the luminosity 
function of galaxies. 

ACKNOWLEDGMENTS 

We thank J.-S. Huang and collaborators for providing 
their K-band luminosity function data in electronic form 
and Simon White and Peder Norberg for helpful conver- 
sations. AJB acknowledges the hospitality of the Institute 
of Computational Cosmology of the University of Durham 
where much of this work was completed. RGB acknowl- 
edges support from a Leverhulme Fellowship. This work 
was supported in part by the UK Particle Physics and 
Astronomy Research Council, the Royal Society and the 
EC Research Network for Research into the Physics of the 
Intcrgalactic Medium. 



REFERENCES 



Adclberger K. L., Stcidcl C. C, Shaplcy A. E., Pettini M., 2003, 

ApJ, 584, 45 
Allen S. W., et al., 2001, MNRAS, 324, 842 

Allen S. W., Schmidt R. W., Fabian A. C, Ebeling H., 2003, 

MNRAS, 342, 287 
Andreon S., 2002, A&A, 382, 495 
Baugh C. M., 1996, MNRAS, 280, 267 

Benson A. J., Cole S., Frenk C. S., Baugh C. M., Lacey C. G., 2000, 

MNRAS, 311, 793 
Benson A. J., Pearce F. R., Frenk C. S., Baugh C. M., Jenkins A., 

2001, MNRAS, 320, 261 
Benson A. J., Lacey C. G., Baugh C. M., Cole S., Frenk C. S., 2002, 

MNRAS, 333, 156 
Bcrlind A. A. et al., 2003, ApJ in press (astro-ph/0212357) 
Blanton M. et al., 2003, ApJ, 592, 819 

Bohringer H., Matsushita K., Churazov E., Ikebe Y., Chen Y., 2002, 
A&A, 382, 804 

Bacon D., Massey R., Refrigicr A., Ellis R., 2002, astro-ph/0203134 
Bower R. G., Benson A. J., Lacey C. G., Baugh C. M., Cole S., 

Frenk C. S., 2001, MNRAS, 325, 497 
Briiggen M., Kaiser C. R., Churazov E., Enfilin T. A., 2002, MNRAS, 

331, 545 

Cavaliere A., Lapi A., Menci N., 2002, ApJ, 581, LI 
Cole S., 1991, ApJ, 367, 45 



Cole S., Aragon-Salamanca A., Frenk C. S., Navarro J. F., Zcpf S. E., 

1994, MNRAS, 271, 781 
Cole S., Lacey C. G., Baugh C. M., Frenk C. S., 2000, MNRAS, 319, 

168 

Cole S. et al., 2001, MNRAS, 326, 255 
Cooray A., Sheth R., 2002, Ph.R., 372, 1 
Cowie L. L., McKee C. F., 1977, ApJ, 211, 135 
Efstathiou, G., 2000, MNRAS, 317, 697 
Eke, V.R., Navarro, J.F. & Frenk, C.S. 1998, ApJ, 503, 569 
Eke, V.R., Cole, S.M. & Frenk, C.S. 1996, MNRAS, 282, 263 
Enfilin T. A., Wang, Y., Nath B., Biermann P. L., 1998, A&A, 333, 
47 

Ettori S., De Grandi S., Molendi S., 2002, A&A, 391, 841 

Fabian A. C, Mushotzky R. F., Nulsen P. E. J., Peterson J. R., 2001, 

MNRAS, 321, L10 
Fabian A. C, Voight L. M., Morris R. G., 2002, MNRAS, 335, 71 
Gruzinov A., 2002, astro-ph/0203031 

Helly J. C, Cole S., Frenk C. S., Baugh C. M., Benson A. J., 

Lacey C, Pearce F. R., 2003, MNRAS, 338, 913 
Hockstra H., Yee H.K.C., Gladders M.D., 2002, ApJ, 577, 595 
Huang J. S., Glazcbrook K., Cowie L. L., Tinney C., 2002, astro- 

ph/0209440 

Jarvis M., Bernstein G. M., Fischer P., Smith D., Jain B., 
Tyson J. A., Wittman D., 2003, AJ, 125, 1014 



12 



Benson et al. 



Jenkins A., Frenk C. S., White S. D. M., Colberg J. M., Cole S., 
Evrard A. E., Couchman H. M. P., Yoshida N., 2001, MNRAS, 
321 372 

Johnstone R. M., Allen S. W., Fabian A. C, Sanders J. S., 2002, 
MNRAS in press 

Murali C, Katz N., Hcrnquist L., Weinberg D. H., Dave R., 2002, 
ApJ, 571, 1 

Kauffmann G., White S. D. M., Guiderdoni B., 1993, MNRAS, 264, 
201 

Kauffmann G., Colberg J. M., Diaferio A., White S. D. M., 1999, 

MNRAS, 303, 188 
Kay S. T., Pearce F. R., Frenk C. S., Jenkins A., 2002, MNRAS, 

330, 113 

Kochanek C. S., Pahre M. A., Falco E. E., Huchra J. P., Mader J., 
Jarrett T. H., Chester T., Cutri R., Schneider S. E., 2001, ApJ, 
560, 566 

Kennicutt, R. C, 1983, ApJ, 272, 54 

Lacey C, Guiderdoni B., Rocca-Volmerange B., & Silk J. (1993), 

ApJ, 402, 15 
McNamara B. R. et al., 2001, ApJ, 562, 149 
Martin C. L., 1999, ApJ, 513, 156 

Mathewson D. S., Ford V. L., Buchhorn M., 1992, ApJS, 81, 413 
Mclchiori et al., 2002, astro-ph/0212276 

Navarro, J.F., Frenk, C.S. & White, S.D.M., 1996, ApJ, 462, 563 
Navarro, J.F., Frenk, C.S. & White, S.D.M., 1997, ApJ, 490, 493 
Narayan R., Mcdvcdev M. V., 2001, SpJ, 562, L129 348, 711 
Norberg P. et al., 2002, MNRAS, 336, 907 

Nulscn P. E. J., David L. P., McNamara B. R., Jones C, 

Forman W. R., Wise M., 2002, ApJ, 568, 163 
O'Meara J. M., Tytler D., Kirkman D., Suzuki N., Prochaska J. X., 

Lubin D., Wolfe A. M., 2001, ApJ, 552, 718 
Peacock, J. & Smith, R. 2001, MNRAS, 318, 1144 
Peacock J. A., 2003, AIP Conf. Proc. 666, 275 



Pearce F. R., Jenkins A., Frenk C. S., White S. D. M., Thomas P. A., 
Couchman H. M. P., Peacock J. A., Efstathiou G.„ 2001, MNRAS, 
326, 649 

Peterson J. R. et al., 2001, A&A, 365, L104 

Pcttini M., Rix S. A., Steidel C. C, Adelberger K. L., Hunt M. P., 

Shapley A. E., 2002, ApJ, 569, 742 
Quilis V., Bower R. G., Balogh M. L., 2001, MNRAS, 328, 1091 
Rees M. J., Ostriker J. P., 1977, MNRAS, 179, 541 
Shu C, Mo H. J., Mao S., astro-ph/0301035 
Sievcrs J. L. et al., 2003, ApJ, 591, 599 
Smith G., et al., 2003, ApJ, submitted. 
Somervillc R. S., Primack J. R., 1999, MNRAS, 310, 1087 
Spergel, D. et al. 2003, submitted to ApJ, astro-ph/0302209 
Spitzcr L., 1962, Physics of Fully Ionized Gases, New York: Wiley- 

Intcrscicucc 
Springel V., Hernquist L., MNRAS, 339, 289 
Strickland D. K., Stevens I. R., 2000, MNRAS, 314, 511 
Tamura T. et al., 2001, A&A, 365, L87 

Taylor G. B., Fabian A. C, Allen S. W., 2002, MNRAS, 334, 769 
van Kampen, E., Jimenez, R., Peacock, J. A., 1999, MNRAS, 310, 
43 

Vianna P.T.P., Nichol R.C., Liddle A.R., 2002, ApJ, 569, L75 
Voit G. M., Bryan G. L., Balogh M. L., Bower R. G., 2002, ApJ, 
accepted 

Voigt L. M., Schmidt R. W., Fabian A. C, Allen S. W., 

Johnstone R. M., 2002, MNRAS, in press 
Wang X., Tegmark M., Jain B., Zaldarriaga M., 2003, astro- 

ph/0212417 

Weinberg D. H., Dave R., Katz N., & Hernquist L., 2003, submitted 

to ApJ, astro-ph/022417 
White S. D. M., Rees M. J., 1978, MNRAS, 183, 341 
White S. D. M., Frenk C. S., 1991, ApJ, 379, 52 
Yoshida N., Stoehr F., Springel V., White S. D. M., 2002, MNRAS j 

335, 762 



APPENDIX 
EJECTION AND RECAPTURE OF SUPERWINDS 

The wind is ejected at speed Vsw, so has mean energy per unit mass of E av = i V S 2 W . We assume that to escape the halo it 
needs energy per unit mass Asw^isk- Assuming a thermal distribution of energies in the wind (i.e. f(E) oc exp(— E/E av )), 
the fraction of the mass flowing out of the disk which also escapes from the halo is given by 

poo 

/eject = / cxp(-.T)da; = cxp(-a;sw), (Al) 
J isw 

where x = E/E av and xsw = AswKnsk/^'av- For the material which escapes the halo, the mean energy excess over that 
needed to eject it from the halo is 

/>00 I POO 

E csc — E av / (x - x S w) exp(~x)dx / / exp(-x)dx — E av , (A2) 

i.e. the mean (kinetic+thermal) energy per unit mass for the gas escaping from the halo is the same as the mean 
(kinctic+thermal) energy per unit mass for the gas flowing out of the disk (since only the highest energy particles escape). 
Using therefore the same energy distribution, we assume that the fraction of particles recaptured in a halo of virial velocity 
^haio is 

/cap = / cxp(-ir)dx = 1 - exp(-a; cap ) (A3) 
Jo 

where x cap = C ca pV^ alo /E av , and C cap is a numerical coefficient of order unity. We have used C cap = 1. 



